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Abstract. 

The inclusive production of several jets is studied in the framework of the hard 
pomeron model. The average jet momenta are found to be strongly ordered and growing 
towards the central cm. rapidity, as expected. Strong positive correlations are found 
for pairs of jets neighbouring in rapidity, both in the forward-backward and forward- 
forward directions. 
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1 Introduction. 

Recent experimental data on low x behaviour of the proton structure functions at 
HERA [ 1 ] seem to support the idea that this behaviour is governed by an exchange of 
a hard (" BFKL" ) pomeron. However a completely inclusive character of the structure 
function measurements makes it very difficult to draw definite conclusions on this point. 
More information is likely to be obtained from less inclusive experiments in which the 
structure of the produced particles spectrum predicted by the hard pomeron model 
could be observed. In a recent publication the exclusive production rate for n jets 
has been calculated in this model [ 2 ]. Since it is not infrared stable, a cutoff on the 
lower limit of the transverse jet momentum had to be introduced. As a result, analytic 
calculations became impossible. Numerical results presented in [ 2 ] seem to depend 
crucially on the introduced infrared cutoffs both for real and virtual gluons. 

In this note we draw attention to the fact that the n-jet inclusive cross-sections in 
the hard pomeron model admit analytic study. We calculate these cross-sections with- 
out any infrared regularization and find them infrared stable, as expected. Of course, 
in the infrared region the model cannot pretend to describe realistic QCD phenomena. 
So, to apply the model to the latter, one has to exclude this region in some way, which 
is important for integrated quantities, such as rapidity distributions or multiplicity mo- 
ments. To this purpose introduction of some sort of cutoff is unavoidable. We however 
find that the model itself provides for a certain infrared regularization (although its 
origin is not quite conclusive, as will be discussed in the following), so that it leads to 
finite rapidity distributions and multiplicity moments. 

Commenting on the results, we stress a characteristic dependence of the n-jet pro- 
duction rate on the jet transverse momenta kf. Since the model is essentially logarith- 
mic, it predicts that all Q = In k\ have the order of at least y/Y where Y = In s is the 
overall rapidity interval. It also predicts that all the differences Q — Q have the same 
large order of magnitude. This means that in the model jet configurations are favoured 
with transverse momenta of a widely different magnitude. Due to dispersion in Ink 2 , 
the momenta of the observed jets result strongly ordered, growing from the projec- 
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tile or target regions towards the center region. Also strong positive correlations are 
predicted both in the forward-backward and forward-forward directions. This features 
could be used to seek for an experimental evidence of hard pomerons. 

The contents of this note is arranged as follows. In Section 2 we deduce a general 
formula for the n-jet inclusive cross-section as a function of jet rapidities and transverse 
momenta in the form suitable for further analytic study. In Section 3 we study the single 
and double jet inclusive cross-sections in some detail, having in mind their practical 
importance. Section 4 is dedicated to a discussion of higher order inclusive cross- 
sections, mostly in certain asymptotic regions of the phase space. Section 5 presents 
some conclusions. 



2 The n-jet inclusive cross-sections 

In the framework of the hard pomeron model jets are currently identified with emitted 
real gluons. The n-jet inclusive cross-section is then described by the diagram shown 
in Fig. 1, where also some notations for the transverse momenta and their relative 
energetic variables are introduced. The blobs marked G denote pomeron forward Green 
functions which describe its propagation between neighbouring jets and also between 
the jets with a maximal or minimal rapidities and the projectile or target, respectively. 
The latter are characterized by their respective colour densities p\ and p 2 . 

Emission of a real gluon (jet) of transverse momentum k is described by a vertex 

l2as -l 2 l f2 (2n) 2 5 2 (l-l'-k) (1) 



k 2 

where / and I' are the reggeized gluon momenta before and after emission; a s is the 
(fixed) strong coupling constant. Passing to the coordinate space, this vertex trans- 
forms into 

^ A «* A (2) 

where the Laplacian operators act in the directions indicated by the arrows. Combining 
these vertices with the pomeron Green functions and densities p we obtain for the n-jet 
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inclusive cross-section 

G„(r,n) A t e*" A, G si2 (r 1 ,r 2 ) A 2 ... A„ G s „ +1 (r„.r') (3) 

Here s^+i denotes the energetic variable for the pomeron between the i-th and % + 1-th 
jets. In terms of jet rapidities yi and transverse momenta fcj 

s iji+ i = I^H/ci+ie^"^ 1 , i^0,n, 

soi = v / s|^i|e~ !/1 , s n , n+ i = v / s|fcn|e I/n (4) 

Evidently they satisfy 

n n 

n = s n *i ( 5 ) 

i=0 i=l 

Eq. (3) is valid provided all energetic variables between jets are large and all momenta 
squared are much smaller: kf, kf +l « s iji+ i — > oo . Then we approximately have 

In s M+ i =yi- y i+ i = y i>i+ i (6) 

for all ? = 0, ..n, with yo = ~Un+i = In = (l/2)y. Note however, that exponentiat- 
ing (6) reproduces only the dominant factor s on the right-hand side of (5). 

The BFKL pomeron forward Green function G s (r,r f ) is well-known [ 3 ]. In the 
limit s — > oo we can retain only its dominant, isotropic, part 

G (r r') = — [ (Ly*» (7) 

s[ ' } 32tt 2 7 (z/ 2 + l/4) 2 V J 1 ' 

Alternatively we may restrict ourselves to studying only the isotropic part of the n-jet 

cross-section obtained after the integration over all azimuthal angles of fcj, which results 

in taking only the isotropic part (7) of all Green functions. Applying the Laplacian 

operators we find 

1 r dv^^ r 
and a similar formula for the Laplacian operator applied to r'. As a result 
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Putting this into (3) we additionally introduce integrations over n + 1 variables i^j+i, i = 
0, ...n. 

At each emission vertex we then find an integral over the transverse distance r 
between the gluons 

f £ -2-2U, £ ikr = _J[j_( fc/2) 2i, r(l-^) ( } 

J u + iO K 1 J r(l + iu) K ' 

where for the i-th jet v = Vi_x,i — The n-jet cross-section takes the form 

T , ,2s A d 3 4 2 " 12ia, 

) = H T-7l2 a = 7T a * 11 



= i dyidkf 9 4 ^ 8tt 2 ^ 
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/n 
d 2 rd 2 r'p 1 {r)p 2 {r')rr' f] dv i ^ 1 s v ^ — 
i=0 v 1 



l/2 + iu m f(l/2-iu^ n+1 ) 
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Now we recall that the expression (11) is valid only for large enough relative energies 
s i>i+ i. The exponent ui{y) has a maximum at v — 0: 

= A - az/ 2 + C(z/ 4 ) (12) 

where 

A = (12a s /n) In 2 (13) 

is the pomeron intercept and 

a = (42a s /7r)C(3) (14) 

Evidently at large s^j+i small values of all i^i+i give the dominant contribution. Con- 
sequently we can neglect i/s in all places where it is not accompanied by some large 
factors. The two integrations over r and r' then give average transverse dimensions of 
the projectile and target 

n>i(02) = J d 2 rrp 1(2) (r) (15) 

Representing 

. r°° da Jau 



v + iO Jin k 2 2tt 



i / — e 
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we can integrate over all z/'s introducing n integrations over ctj, % = 1, ...n instead. 
Then we obtain for the n-jet cross-section 

y i=l 47r K « i=l i=0 

where 

and we have used the notation 

Ci = h*. ? (18) 

and have assumed ao = a n+ ± = 

Expression (16)-(17) is still far from being explicit, since integrations over a's cannot 
be easily done, except for the simplest case n = 1. However the integrand is now a 
positive function well behaved at all values of cc's. Therefore Eqs. (16)-(17) can be 
conveniently used both for numerical calculations and analytic study in certain limiting 
regions of jet parameters. Before we pass to the latter study we want to point out some 
evident general features of the n-jet cross-section (16)-(17). 

It is clear that its overall behaviour in s is independent of n and is the same as 
for the total cross- sect ion, namely, ~ s A /Vhi s. Indeed the characteristic values of 
all a's are of the order ^JylJ^[ ~ Vhi s. Therefore integrations over a's will give a 
factor ~ (Ins)™/ 2 , which will cancel a similar factor in (16) leaving one Vln s in the 
denominator. Passing to the behaviour in k%, that is, in Q, we observe that at large Q 
the exponentials in the integrand provide for a damping factor at Q > \/Y . As a result 
all C's and also their differences \Q — (j\ have large average values of the order \/Y . This 
leads to a distribution in kf of the form discussed in the Introduction: neighbouring 
jets prefer to have their transverse momenta widely different in magnitude. 

Of special interest is the infrared behaviour of (16)-(17), when all Q — > — oo. The 
integrals over a's converge to a finite value in this limit. Therefore the behaviour of 
the cross-section in each kf will be determined by the factor kf A /kf. Taken at its face 
value, it leads to the cross-section (16) integrable at small kf and thus infrared stable. 
However we have to acknowledge that retaining the factor kf A in our formulas is in 
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fact beyond the accuracy adopted in the lowest order hard pomeron model, in which 
terms containing powers of a s In k 2 have been systematically neglected. Actually the 
region of small kf is not only beyond the applicability of the hard pomeron model to 
the realistic QCD but also beyond the validity of the model itself, since, as follows 
from (4), it corresponds to relative energies s^+i becoming small, when the lowest 
order approximation fails. For this reason the small kf region has to be excluded in 
any case, whether one wants to apply the model to the realistic QCD or to study the 
model as it stands. The factors k A which have appeared in (16) serve precisely to this 
purpose: they effectively cut off the small kf region. Of course, one may alternatively 
use harder cutoffs, as currently applied in the hard pomeron calculations, when the 
integration over k 2 is simply cut off at some minimal fc^ ira . Since the difference, if any, 
reflects only the unknown influence of the infrared region, in the following we shall 
retain the factors k 2A in (16) as a natural infrared cutoff. 
Finally we note that J n , Eq.(17), can be rewritten as 

n -.qc n n—1 

Jn = Y[ci da>iexp(-J2a 2 - 2 b j a j a j+i) ( 19 ) 
i=i Hi/ci j = i j = i 

where 



V yi-i,i+i ' 1 \j yi-i,i+iVi,i+2 

and we have denoted y^i+z = yi,i+i + Vi+i,i+2- From this expression it is clearly seen 
that correlations are predicted between jets neighbouring in rapidity (i.e. i and i + 
1). Their strength is determined by the value of 6-,'s, which depends on the rapidity 
distances between the correlating jets and their neighbours. Evidently the correlations 
are maximal if the jets are emitted at the center of the rapidity interval between 
their neighbours and vanish if they are emitted at its borders. One also sees that the 
correlations are positive if the two jets have large momenta (Q, Q + i > 0), whereas they 
may be both positive and negative if one or both jets have very small momenta (Q 
or/and Q +1 < 0). Study of the 2-jet correlations shows that in the latter case the 
correlations are negligible. 
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3 Single and double inclusive cross-sections 

Although single jet inclusive cross-sections have already been studied rather extensively 
since long ago, we found it convenient to start with them to clearly see the effect of 
the natural cutoff provided by the k 2A factor in (16). With it the production rate can 
be written explicitly. Earlier studies mostly used harder cutoffs, like just excluding the 
region below some chosen minimal value &^ n . 

The single jet inclusive cross section, obtained from (16) for n — 1 can be written 
in terms of the error function 

/i(l/i,*?) = ^fc? A - 1) (l-$(x 1 )) (20) 
where 



Xl V 4c "/0l2/l2 ^ 

and a is the total cross section corresponding to a pomeron exchange, which is a 
common factor to all inclusive cross-sections: 

For finite k\ and thus finite £i, the argument x\ is small and the behavior of I\ in k 2 is 
totally determined by the factor 1/kf. At very small kf (d — > — oo) the cross-section 
is cut by the factor k\ A at values Ci ~ — 1/A. At very large k\ a cutoff is provided by 



the error function at (i ~ ^ 4ayoiZ/i2 / 'Y 

Integrating (20) over kf one obtains the distribution in rapidity 

h(yi) = ^^aexp(aA 2 y 01 y 12 /Y) (23) 

Since yoi2/i2 = Y 2 {\ — y\) this distribution is evidently peaked at the center rapidity. 

One should note that the exponent in (23) has the order a 3 s Y and thus, strictly 
speaking should be neglected in the lowest order approximation, in which only terms 
of the order a s Y has been taken into account. Keeping the exponent has the same 
meaning as retaining the factors k 2A in the initial formulas and should be considered 
as a natural cutoff provided by the kinematical constraints in the theory. 
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To find an average value of the jet momentum at a given rapidity one has to integrate 
(20) with the weight Ci and divide by (23). One then obtains 

(d) = (l/2)oA(y - Ayl/Y) (24) 

As expected from the BFKL diffusion, the average momentum grows from the edges of 
the rapidity interval towards the center. The characteristic order of (i = In k\ results 
to be a 2 s Y, and not \/Y as deduced from the BFKL wave function properties. This 
also is a consequence of retaining the k 2A factor, which evidently somewhat enhances 
the observed jet momentum. 

Integrating (23) over all rapidity range one obtains the total multiplicity as 



{n) = 2SW v exp(aA2y/4)$ G Av/ ^) (25) 

For AV aY << 1 the average multiplicity grows linearly with Y, as expected. However 
if one takes Y — > oo keeping A fixed, (n) rises as a power of energy ~ s aA2 / 4 . 
Now we pass to the double inclusive cross-section given by (16) for n = 2: 



h{yi,kl-y 2 ,kl) = ^a(hk 2 ) 2A -\ ^—J 2 (26) 

167T 5 a V 2/012/122/23 



oc 



J 2 =l daW da 2 e W (--^- {a \- a2)2 - (27) 



where 

ki ™ 1 ^2 "^""^ 4 a y 01 4 a2/l2 4gm/ 2 3, 
As well-known, this cross-section was studied in [4] under different kinematical condi- 
tions (small yoi and y 2 % and large y± 2 ). Only the dependence on y i2 was discussed as a 
possible experimental signature of the BFKL pomeron between the two jets. 

As seen from (19) real parameters of the integral (27) are £i/y/4aY and ( 2 /\^AaY. 
So for finite jet momenta one can put (i = ( 2 = and the integral, which turns out 
independent of jet momenta, can be easily calculated; 



■h = -Ka\jy m y 12 y 2 j,jY (^- + ~ arctan 



/ 2/012/23 
2/12 V\ 



so that the inclusive cross-section becomes 



I 2 ( Vl , h\- y 2 , hi) = M^J^a-i (l + l arctan (28) 
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For general £i and ( 2 the integral (27), which is an evident generalization of the 
error function to two dimensions, cannot be done analytically and requires numerical 
calculation. We have calculated the ratio 



0^2(2/1, fcii 2/2, fcj) 



fci; 2/2, fc 2 ) = T77. ,„2\t 77. 17T\ ( 29 ) 



to study correlations between the emitted jets. Absence of correlations corresponds to 
R = 1. In Fig. 2 we show results of the numerical calculation of R for a symmetric 
pair of jets at y\ = —y 2 = Y/6, which divide the rapidity interval into three equal 
parts. R is presented as a function of z 2 = C2/ V^aY in the interval < z 2 < 2 for fixed 



positive values of z 1 = d/V4aY < 1. For negative Ci or/and ( 2 we find that R = 1 
and there are no correlations. One observes from Fig. 2 that there are significant 
positive correlations between the two jets in the whole region of positive Ci and ( 2 
corresponding to momenta which are not small relative to the intrinsic scale of the 
BFKL model (which should be supplied by higher orders of the theory). They should 
be seen experimentally as positive forward-backward correlations in the cm. system. 

For very large jet momenta, when £1 or/and ( 2 become of the order V 'AaY , the 
correlations become very strong and grow very fast with £'s, reaching values of the order 
4. 10 5 at Ci = C2 = 2. However one should have in mind that at such high momenta 
the production rate gets exceedingly small, beyond any possibility of experimental 
detection, which can be deduced from (20), so that such high values of R bear, in fact, 
no physical importance. 

We have also calculated the ratio R for two jets emitted both in the forward hemi- 
sphere, with y\ = Y/3 and y 2 = Y/6. The results are quite similar qualitatively, 
although the numerical values result a bit different (however of the same order of mag- 
nitude). Thus one should also expect forward-forward correlations of the same order 
as the forward-backward ones. 

Integrating (26) over the two momenta k\ and k\ one obtains the double distribution 
in rapidity. These integrations can easily be done explicitly using an evident formula 

r- + 00 

r/n/(n)= (1/A) / 



/+00 r+00 r+00 

d(e AC / daf(a) = (1/A) / dCe A <f(C) (30) 
-00 JC J -co 
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valid for any f(a) integrable at a = ±00. We obtain 

HviiVz) = (2^;) aexp ( aA2 (yo^ + 2/23 - (2/01 - V23) 2 /Y)) (31) 

The exponential factor again comes from the kinematical factor k 2A and has a 
higher order than normally taken into account in the BFKL theory. If we neglect it, 
we obtain a completely flat distribution in the two rapidities. The exponential factor 
evidently has its maximum when the two jets are both emitted at central rapidity 
Vi — U2 = 0, when it grows like exp(aA 2 V). Comparing to (23) we observe that it also 
introduces positive correlations between the two jets in the rapidity space. In fact, the 
ratio analogous to (29) results 

R(Vi, V2) = f ff;f ] , = exp(aA V W Y) (32) 

h{yi)h{y2) 

and evidently implies positive correlations, growing as the two jets move to the central 
rapidity region. 

Integrating (31) over all the rapidity region we obtain the second moment of the 
distribution in the number of jets. Neglecting the exponential factor we get 

<n(n-l)>=(|^)V = <n> 2 (33) 

With the exponential factor taken into account, the second moment can be found only 
by numerical integration. 

4 Multijet inclusive cross-sections 

The features found for the production of two jets repeat themselves in the multijet 
production. However for n > 2 analytic results cannot be obtained for most of the 
observables. So one can only study cases corresponding to some limiting parts of the 
phase space. 
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4.1 Distributions in rapidities and momenta 

As mentioned the real parameters of the cross-section are ratios d/^/AaY. Depending 
on their values we may distinguish between three asymptotic regions, in which the 
cross-section (16) can be estimated by standard asymptotic methods. 

The first region is the infrared one, when Q — > — oo faster than VY — > oo. In this 
region the cross-section is given by (16) with the integrals J n independent of kf and 
given by 

J n = ft r da, exp (- ± ^- a ^A (34) 
t=i J -°° V j^o 4 «%j+i J 

To do this integral we pass to variables 

a iji+ i = OLi - a i+1 , i = 0,...n (35) 
subject to condition J27=o a i,i+i — an d present 

/oo A \ n roc n ( ty? ■ - \ 

d^iexpE ( iX -Jzf^) (36) 

All integrals are trivial and we get 

^ = ^^n(4vra 2/M+1 ) 1/2 (37) 
Thus in the infrared region we have the cross-section 



,3a. 



n 



Jn(w,A?) = (2^)MlI*?r" 1 (38) 
i=i 

As one observes, the infrared cross-section is independent of individual jet rapidities: 
the distribution in all is completely flat. The dependence on kf apart from the 
canonical factor 1/kf exhibits an extra factor discussed above kf A , which serves as an 
effective infrared cutoff. 

The second region is that of low kf, when Q are finite and thus all ratios C/VY go 
to zero. In this region the inclusive cross-section is given by (16) with the integral J n 
again independent of kf : 

\2 ' 



J 1 roo I n 

i=l J0 \ j=0 



o 4ay jd+1 



(39) 
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The dependence on kf is thus the same as in the infrared region. However there appears 
now a nontrivial dependence on rapidities. 

The integral (39) can be easily done only for n = 1,2 (see Sec. 3). Calculation of 
J n for larger n can be done by numerical methods. Analytic estimate of J n is possible 
within the saddle point method, presumably valid for large n. Denote 

n I \2 n 

p = T, 1 \ = E - «m) 2 ( 4 o) 

i=0 ^ a Ui,i+l i=0 

We have 

fc+m 

E dP/dai = 2b kyk+1 (a k+1 - a k ) + 2b k+m)k+m+1 (a k+m - a k+m+1 ) (41) 

i=fc+l 

If we require that 

dP/dat = 0, (42) 

for all i = 1, ...n then, recalling that ao = a n+ \ = 0, we shall obtain an homogeneous 
system of linear equations in a^i = l,...n with a nonsingular matrix, so that the 
solution will be that all a are zero, i.e. at the initial integration point. Thus there 
are no saddle points inside the overall integration volume in all en's. We can, however, 
relax condition (42), by fixing one of the cc's, say a n , thus seeking a saddle point on 
the surface, inside the integration volume in all other ctj, % = l,...n — 1. Putting 
k + m = n — 1 in (41), we get 

a k - a k+1 = _ ^"~ 1 ' n (an-! - a n ) (43) 

Summing these equations for k — 0, ..j — 1 we obtain 

i-i j, 

- (Xj = 6 n _i ;n (a; n _i - «„) E V^.fe+i = '" K-i - «n) (44) 

fc=0 °0j 

Here and in the following we use the notation 

j+m-l 

l/b jJ+m = E V 6 fc,fc+1 = 4a (2/j - %+m) ( 45 ) 
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Putting j = n in (44), we can express the difference a n -i — a n in terms of the fixed 
variable a n . Combining this with (44) we finally find 

<Xj = - — a n (46) 
b o,j 

Eqs. (46) determine the position of the saddle point in a iy i — 1, ...n — 1 in terms of 
the fixed a n . At the saddle point we find 

n-l _ \ 2 

P =J2 Ki+i H~ 1,n ( a n-i - Oi n ) + b n<n+1 a 2 n = a 2 n (b 0n + 6 n , n+ i). (47) 
i=o V / 

Integration over the saddle point gives a factor 



n n-l 



V detM 

where M is an (n — 1) x (n — 1) matrix of second derivatives. Simple calculations give 

detM = (l/6o,„) II Ki+i ( 48 ) 

i=0, n— 1 

The integration over a n adds a factor lyV/ (6 0n + ). and we get the contribution 
from the surface a n = 



1 / Tjf. n 

^=^V^n(4^ + i) V2 (49) 

which is one half of (37). Taking into account all other surface contributions will 
multiply (49) by n. As a result, we obtain for the inclusive cross-section an expression 
which is identical to (38) except for an extra factor n/2. The distribution in rapidities 
obtained by this asymptotic estimate is again completely flat. 

We finally pass to the third asymptotic region, that of high kf, which corresponds 
to all Q/ y/y — > oo. In this region the behaviour of the integral (17) is governed by the 
damping exponentials in the integrand. To find the asymptotics we may again recur 
to the saddle point method. However in this region the integration volume depends on 
the values of Q and we have to check whether the found saddle point lies inside this 
volume. As a result, the set of a's for which the saddle point lies inside the integration 
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volume will depend on the values of Q and will generally be only a part of all cc's. The 
rest of a's will be restricted to their initial values £'s, since in absence of a saddle point 
the integrand falls with each a exponentially. The asymptotics will thus depend in a 
rather complicated way on the manner in which different £'s get large and also on the 
jet rapidities. Nevertheless it is not difficult to find this asymptotics for given yi and 
large Q. 

Let us consider a situation when a particular set of a's take on their initial values 
in the asymptotics, the rest of them integrated over their saddle points. Take a pair 
of jets, the j th and (j + m)th, for which cc's take on their initial values with all 
intermediate jets from (j + l)th to (j + m — l)th having a saddle point inside the 
integration volume. In other words Eqs. (42) are satisfied for i — j + 1, ...j + m — 1 
with otj = Q and aj +m = Cj+m- Repeating the procedure described above we can solve 
these equations for the saddle point ctj, i — j + 1, ...j + m — 1 in terms of Q and (j+m-i- 
We find 

a j+ k = Q + -r^-iCj+m ~ Cj) > Cj+k, fc=l,...m-l (50) 

Oj,j+k 

The last inequality follows from the condition that the saddle point should lie inside 
the integration volume. This inequality determines a certain domain in the jet phase 
space i/i, kf where this particular type of asymptotic behaviour will take place. The 
part of the exponential P containing a's from the (j + l)th to (j+m — l)th, at the 
saddle point becomes equal to 

that is, takes the same form as if there were no other jets between the j th and (j + m) 
th ones. Integration over all aj, i — j + 1, ...j + m — 1 by the saddle point method thus 
gives a factor 



^m—lU. . 
\ lli=j 

and leaves an integrand where these jets have disappeared. Doing this with all sets of 
a's which have their saddle point inside the integration volume we shall obtain more 
factors (51) and an integral over remaining a's in which all other jets have disappeared 
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and which does not contain saddle points inside the integration volume. Integrating 
by parts one obtains for the latter integral an asymptotics 



Both the product and sum go over the jets remaining after exclusion of the ones inte- 
grated over their saddle points. The described asymptotics will evidently be valid in 
the domain where the inequality (50) is fulfilled for these excluded jets. 

To make this rather complicated prescription more understandable, consider a sim- 
ple case when rapidity distances between neighbouring jets are all equal: y^i+i = Ay 
independent of i. Then the left hand side of inequality (50) represents a direct line 
in the i, ( plane connecting points Q and Q +m . The inequality (50) means that the 
intermediate Q, % — j + 1, ...j + m — 1 lie below this line. The prescription for the 
asymptotics is then to draw all lines connecting given Q, i — 0,1, ..n, n + 1, where by 
definition Co = C«+i — 0- If some £'s happen to lie below any such line, one should 
throw the corresponding jets out introducing factors (51) instead. In the end only C's 
will remain which lie above any line connecting a pair of them, that is the remaining 
Q will represent a convex line in the i, ( plane. The behaviour in these C's will be de- 
scribed by (52). Of course, if if all Q form a convex line from the start, the asymptotics 
will have the form (52) with a product and sum going over alH = 1, ...n. 

For illustration consider a case n = 2. Condition that points Co — 0,Ci,C2,C3 = 
form a convex function of % is evidently 




(52) 



ICi - C 2 



< Ci, c 2 



(53) 



In this region the asymptotics will be 



1 



-exp 



( 



Ci 2 + CI + (Ci 




) 



(54) 



47T 5 



C1C2 " 2(Ci - C2) 



(4/3)aF 



Now take the region 



2Ci < C 2 



(55) 
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where (j lies below the line connecting (o = and £2- According to our prescription 
the first jet now gives a factor (51) and the second gives (52). We thus obtain 



2 



^^-V^-'^pf-^j (56) 

In the symmetric region 2( 2 < Ci we shall evidently obtain (56) with the first and 
second jets interchanged. 

4.2 Integrated distributions 

Integrating (16) over all kf we get rapidity distributions I n (yi). As discussed, these 
integrations converge in the infrared region due to the cutoff factors kf A . Using (30) 
we obtain from (16) 

UVi) = ^n^T-VYlly^^ (57) 

where 

/,+oo n / n n \ 

IlrfCexp A£0-EW0-0 + i) 2 (58) 
-°° i=l \ i=l i=o / 

and Co = Cn+i — an d h,i+i are defined by (40). The integral (58) can be easily 
calculated using the same technique that was applied to obtain (37). In terms of 
variables C»,*+i = d ~ Ci+i> the sum °f which is zero, we find 

n n 

£Ci = (l/2)£(n-2i)C i)i+ i (59) 

i=l i=0 

The integral K n becomes 

/+00 d\ n . . 

— n^OexpE f(«A + (ra/2 - i)A)C M +i - &M+i(Ci " Cm) = 

1 n /aA 2 \ 

— J[(A7ray ht+1 ) l/2 exp (^(n 2 F 2 + 8pY - 4(n + 1)<tY - 4. 2 )j (60) 



where 



fc=i fc=i 
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Putting (60) into (57) we obtain finally 



In{Vi) = ( 




(n 2 Y 2 + 8pY - A(n + l)aY - Aa 2 ) 



(62) 



As in the two-jet case, the exponential factor lies, in fact, outside the precision, 
normally adopted in the lowest order BFKL model. If we neglect it we evidently 
obtain a completely flat distribution in all rapidities. The effect of the exponential 
factor can easily be seen. Let = n 2 Y 2 + 8pY — A(n + \)oY — Aa 2 . Since the 
derivatives d(f)/dyk = Y(8k — A(n + 1)) — 8a, k = 1, ...n cannot all be equal to zero 
simultaneously, the maximum of <fi occurs at limiting values of its arguments, that 
is, when all yi coincide: yi — y 2 — ... = y n . At this point = n 2 (Y 2 — Ay 2 ). The 
maximum is achieved at y\ — and is equal to n 2 Y 2 . Thus the exponential factor 
favours emission of all jets at central rapidity transforming a flat plateau into a convex 
one. 

Integration of I{yi) over all yi gives multiplicity moments. With a flat distribution 
we immediately get 



This corresponds to a Poissonian distribution in n, with the average (n) growing like 
Y, in confirmation of simple estimates. The exponential factor in (62) will, in all 
probability, make the moments grow more rapidly, like powers of energy, as follows 
from the case n — 1 studied in Sec. 3 

5 Conclusions 

We have studied inclusive production of many jets in the hard pomeron model, with 
a special emphasis on the two-jet production. In confirmation of current expectations, 
we found that the average jet momenta are strongly ordered, monotonously growing 
from the projectile or target regions to the central regions. The most characteristic 
feature of multiple jet production seems to be strong positive correlations between jets 
neighbouring in rapidity, both in the forward-backward and forward-forward directions. 



(n(n - l)...(n -& + !)) = ( 



3a s 



) k Y k = (n) k 



(63) 



2tt 2 A 
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In our treatment we have preserved kinematical factors k , which serve as a nat- 
ural infrared cutoff and allow to obtain finite analytic expressions for distributions in 
rapidity and number of jets. As a consequence we also obtained an enhancement in the 
high k 2 region, so that the average values (Ink 2 ) result growing like a 2 Y and not like 
y/a s Y as concluded from the study of the pomeron wave function. To see the relevance 
of this effect, one has to study the second order corrections in the hard pomeron model, 
which now seems possible, but lies beyond the scope of the present study. 
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8 Figure captions 

Fig. 1. The diagram which describes the n-jet inclusive cross-section. 



Fig. 2. The ratio R (Eq. (29)) for two jets with y 1 = —y 2 = Y/3 as a function of 



In hf/V^aY for different values of z\ = In k\/\fAaY . 
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